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Entanglement witness is a Hermitian operator detecting entanglement in a quantum system. Non¬ 
linear entanglement witnesses have wider range entanglement detection ability relative to the linear 
ones. Here we construct some entanglement witnesses for a four qubit mutually unbiased basis Bell 
diagonal density matrices. To this aim we specify the convex feasible region with positive partial 
transpose states. To reveal the entangled regions, we present some appropriate linear entanglement 
witnesses as the envelope of family of linear entanglement witnesses and in this case, the nonlinear 
witnesses have been nonlinear functional of linear witnesses. We also study thermal entanglement 
and show for some Hamiltonians the witnesses have ability to detect the entanglement at all temper¬ 
atures. The results of this study are evaluated and assessed in the light of four qubit entanglement 
detection problem. 

PACS numbers: 03.67.Mn, 03.65.Ud. 


I. INTRODUCTION 

Quantum entanglement is a physical phenomenon that 
occurs when pairs or groups of particles are generated or 
interact in ways such that the quantum state of each par¬ 
ticle cannot be described independently instead, a quan¬ 
tum state may be given for the system as a whole [1, 2]. 
Mathematically, a state of a composite quantum system 
is called entangled if it cannot be written as a convex 
combination of product states Q • Quantum entangle¬ 
ment has many physical application such as quantum key 
distribution in quantum cryptography [4 5] with new ex¬ 
periments @,[3, quantum dense coding [8|, and quantum 
teleportation [§, j39|. In these applications there must be 
some physical observable acting on system state to detect 
the entanglement in the system. On of the observables 
detecting entanglement especially for system with three 
or more particles is entanglement witness (EW). 

Entanglement witness is an observable which com¬ 
pletely characterize separable states and detect entan¬ 
glement in a system experimentally EMI. From ge¬ 
ometrical point of view as the quantum mixed states 
family (density matrices) is a convex set so an EW 
can be describe by hyperplanes in the density ma¬ 
trix space. Now the EW, IT is a Hermitian op¬ 
erator which have a non-negative value on all pure 
product states Tr(Wp Product ) > 0 where p Product = 
| )... | ifri ) {iftn I • • • {Vh I • The entanglement of p is detected 
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by EW if and only if Tr(Wp) < 0 Although there is 
a necessary and sufficient condition for separability in 
2 <g) 2 and 2 <g) 3 cases, called the positive partial trans¬ 
pose PPT criterion or Peres-Horodecki criterion EH, in 
general there is no such conditions for other cases and 
there are states that are entangled but PPT in all those 
cases which are called PPT entangled states. The an¬ 
other way for detecting entanglement for system with 
higher dimensions is using EW. Specially the EWs de¬ 
tecting PPT entangled states are in great importance. 
Usually these are non-decomposable EWs. One can con¬ 
sider linear EWs which is relatively simple to construct 
or nonlinear EWs 

There are some nonlinear separability criteria in the 
literature. Generally these nonlinear EWs have wider 
range entanglement detection. In article Ea , the author 
derive a family of necessary separability criteria for finite 
dimensional systems based on inequalities for variances 
of observables and formulate an equivalent criterion in 
terms of covariance matrices. The criteria may be ap¬ 
plied from the regime of continuous variables to finite 
dimensional systems. Nonlinear EWs as an extension of 
linear witnesses with ability to detect the states with neg¬ 
ative partial transpose has been presented in EH ■ A gen¬ 
eral theorem as a necessary condition for the separability 
based on concave-function uncertainty relations has been 
derived for both finite and infinite dimensional systems 
in EH- The author has been used the specific concave 
function method for a system with mutually unbiased 
bases (MUB), for entanglement detection as a special 
case of his approach. In some cases that approach leads 
to an analytic entanglement detection which is stronger 
than the Shannon entropic uncertainty relation and the 
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Landau-Pollak uncertainty relation. Using appropriate 
class of uncertainty relations, the entanglement for the 
local quantum states of a pair of IV-level systems have 
been defined in [20| . These uncertainty relations may be 
used as an experimental test of entanglement generation. 
A derivation of nonlinear EWs based on covariance ma¬ 
trices has been investigated in [2lj |. The nonlinear func¬ 
tions which improves the entanglement detection given 
by the linear ones is presented in f22j with explicit ex¬ 
amples showing accessible nonlinear EWs detect more 
states than their linear ancestors. 


The other way for constructing nonlinear EWs is based 
on the PPT entangled states detection by improving the 
linear EWs. In this approach for a given density ma¬ 
trix the PPT convex region is determined by the PPT 
inequalities of the density matrix. For some PPT states 
in this region which called feasible region, the PPT cri¬ 
terion is sufficient for separability. Then linear EWs in¬ 
troduce and the non decomposability for some of them 
is achieved. Then nonlinear EWs has been constructed 
from linear ones. This method has been applied for three 
qubit MUB diagonal entangled states in [23j], for 2 ® 2 ® d 
bound entangled density matrices by exact convex opti¬ 
mization in [241 ] . for general algorithm for manipulating 
nonlinear and linear entanglement witnesses by using ex¬ 
act convex optimization [251 ] . and for bipartite N ® N 
systems via exact convex optimization in jfffij. 


In this paper we aim to determine the EWs for a given 
four qubit density matrix which have ability to detect the 
PPT entangled states. Then we are intended to provide 
the nonlinear EWs from the linear ones. To do this first 
we specify the PPT states for a given four qubit density 
matrix in the mutually unbiased basis (MUB). These 
states make a convex region called the feasible region 
( FR ). By using the PPT inequalities we find some MUB 
diagonal states which for them the PPT criterion is nec¬ 
essary and sufficient for separability. Then we introduce 
the linear EWs family which detect the bound MUB di¬ 
agonal density matrices. Then we present the nonlinear 
EWs with the nonlinear coefficients which have wider 
range detection. These nonlinear EWs are envelope the 
family of linear EWs. This framework helps to inves¬ 
tigate the EWs for a given density matrix (here four 
qubit) and serves for in organizing the knowledge about 
entanglement of the system. In the last section we study 
temperature effects on the entanglement for an ensem¬ 
ble of four qubit systems in thermal equilibrium. Finally 
we compare entanglement detection of the nonlinear wit¬ 
nesses with state preparation fidelity. The result shows 
for special case our nonlinear witnesses can detect the 
entanglement at any temperature for some coupling con¬ 
stants while the fidelity method fails. 


II. MUB DIAGONAL DENSITY MATRICES 
AND POSITIVE PARTIAL TRANSPOSE 
CONDITIONS 


Mutually unbiased bases (MUB) in N dimensional 
Hilbert space are orthonormal bases |vi) and | Wj) such 
that |(ui|wj}| = 1/VN for all i, j £ {1, ..,1V}. If one 
can find N + 1 mutually unbiased bases for a complex 
vector space of N dimensions, then the measurements 
corresponding to these bases provide an optimal means 
of determining the density matrix of an ensemble of sys¬ 
tems [27l - |29fl . These bases may be used for entanglement 
detection EkMbil . 

The Bell basis is an orthonormal basis for the two qubit 
Hilbert space and in terms of computational basis could 
be written as 

hr) = -±=(|oi) -| 10 », |V’ + ) = ^(| 01 ) +| 10 », 

i«n = ^=(ioo) - in)), |0+) = -^doo) +in».(i) 

Bell state is any quantum state in the Bell basis. Density 
matrices which are diagonal in this basis are called Bell- 
diagonal. In the case of two qubit the Bell-diagonal state 
is 


=Pl|^ + )(</ ,+ l+P2|V’ + )(V’ + l+P3|V' ){ 4 > I+P4I0 )( 4 > I 

( 2 ) 

where 0 < pi < 1 and Y^i=iPi ~ 1- 

The Bell basis can be generalized, specifically consider 
a system of four qubit spins, the generalized 16 elements 
can be written as 


K) = i=(|oooo)-Hiiii)), |v- 2 > 
|^) = -±=(|0001) + |1110)), |V- 4 > 
| t ) = i=(|ooio)-Hiioi)), |v> 6 ) 
IV> 7 ) = -±=(|0011) + |1100)), |^ 8 ) 
K > = ^=(|0100) + |1011», |^ 10 ) 
|^ 11 ) = -^(|0101) + |1010)), |^ 12 ) 
|V> 13 ) = ^(|0110) + |1001)), |V> 14 ) 
|V> 15 ) = ^(|0111)+|1000», |^ 16 ) 


—(|0000) 11111)) 
-±=(| 0001 >- 11110 )) 
-±=(| 0010 ) — 11101 )) 
-±=(| 0011 )- 11100 )) 
-±=(| 0100 ) — 11011 )) 
-±=(| 0101 )- 11010 )) 
-±=(| 0110 ) — 11001 )) 
-±=(| 0111 )- 11000 )) 


(other basis choices are possible, for example see }33l . l34i ]). 
The diagonal Hamiltonian in this bases is 

16 

h = J 2 e * i^x^i ( 3 ) 

2=1 
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where E z is the energy eigenvalue of the |state. In 
terms of Pauli spin matrices 


H = c 0 IIII + c 1 Ia z a z I + c 2 Ia z Ia z + c 3 IIa z a z + 

c 4 a z IIa z + c 5 a z Ia z I + c 6 cr z a z II + c 7 a z a z a z a z + 
C^O'xO'xO'xO^x d" C 9 @x@y®y@x T Oj^ O' x (JyJ x Jy T 
c^o'xO'xO'yO'y T c i2 a y a x a y a x T c^j y o X JxJ y T 

C i4 <7 y a y IJ x C! 'x + C i 5 a y a y cr y a y (4) 

( the tensor product sign is omitted for simplicity, for 
example the forth term cr z IIa z means a z ® I I ® a z ) 
here c z s can be driven in terms of E z s and characterize 
the coupling strength among qubits. The first term rep¬ 
resent no interaction at all (a constant term), the next six 
terms represent pair z component spin interaction (Ising 
like), and the remaining terms represent the four party 
interactions. 

Now suppose we have a large number (theoretically, 
infinite) of four qubit molecules in thermodynamic equi¬ 
librium (canonical ensemble). If we assume that the in- 
termolecular interactions are negligible, then the total 
system is in a product state, p <8> ■ ■ ■ <8> p, it follows from 
the additive property of entanglement that the total en¬ 
tanglement present in the system is N times the entan¬ 
glement present in a single molecule, where N is the total 
number of molecules present in the system j3f|. 

For a canonical ensemble of four qubits in the thermal 
equilibrium the state of this system in the Bell-diagonal 
bases can be written as 


16 

P = 5> IV’iX^il 

i=1 


here 


Pi = 


0-0 Ei 




( 5 ) 


( 6 ) 


is the probability of finding the system in the state \i/>i), 
and /3 = ks is the Boltzmann constant, T is the 

temperature and Y^]=iPi = 1- 

We can write the density matrix in terms of two di¬ 
mensional Pauli matrices as follows 


p = — [IIII + r 1 Ia z a z I + r 2 Ia z Ia z + r 3 IIa z a z + 
r 4 a z IIa z + r 5 a z Ia z I + r 6 a z a z II + 

T 7 @ Z &Z&Z&Z T T' 3 (7 x JOx J'X T T Q <7x<J y<7 yOx T 
^ 10 U 7 ; Gy G X J/j T l’ 11 C r X C7 a ;(T y(7y “f“ V ^ 2 J y (J x (Jy (7 X T 

T^CTyO'xO'xO'y T 1*! 4 Jy Jy J X J 7 -f- 1", ^ J y Jy Jy Jy ^ (T) 


where the coefficients are 


r l 

= P1 + 

p 2 

+ p 3 

+ P4 

-Po 

— 

Pg 

1 

-I 

00 

— 

P 9 

— 

PlO 


-Pll 

- 

Pl2 + Pl3 

+ Pl4 

+ 

Pig 

+ PlG 






r 2 

= P1 + 

p 2 

~P 3 

-P 4 

+ Po 

+ 

Pg 

-P7- 

Ps 

- 

P 9 

- 

P 10 


+Pll 

+ 

P12 - 

P 13 

P 14 

+ 

Pig 

+ PlG 






r 3 

= P1 + 

p 2 

~P 3 

P 4 

P 5 

- 

Pg 

+ P 7 + 

P 8 

+ 

P 9 

+ 

PlO 


-Pll 

- 

Pl2 - 

Pl3 

P 14 

+ 

Pig 

+ PlG 






r 4 

= P1 + 

p 2 

-P 3 

-P 4 

+ Po 

+ 

Pg 

~P 7 - 

Pg 

+ 

P 9 

+ 

PlO 


-Pll 

- 

Pl2 +Pl3 

+ Pl4 

- 

Pig 

-PlG 






r 5 

= P1 + 

p 2 

+ P 3 

+ P 4 

-Pg 

- 

Pg 

~P 7 - 

Pg 

+ 

P 9 

+ 

PlO 


+Pll 

+ 

Pl2 - 

Pl3 

~Pl4 

- 

Pig 

-PlG 






r G 

= P1 + 

p 2 

+ P 3 

+ P 4 

+ Pg 

+ 

Pg 

+ P 7 + 

Pg 

- 

P 9 

- 

PlO 


-Pll 

- 

Pl2 - 

Pl3 

~Pl4 

- 

Pig 

-PlG 






r 7 

= P1 + 

p 2 

-P 3 

-P 4 

-Pg 

- 

Pg 

+ P 7 + 

Pg 

- 

P 9 

- 

PlO 


+Pll 

+ 

Pl2 + Pl3 

+ Pl4 

- 

Pig 

-PlG 






00 

= Pi - 

p 2 

+ P 3 

-P 4 

+ Pg 

- 

Pg 

+ P 7 - 

Pg 

+ 

P 9 

- 

PlO 


+Pll 

- 

Pl2 +Pl3 

~Pl4 

+ 

Pig 

-Pig 







r 9 = -Pi + P 2 - P 3 + P* + Ps - Pe + P 7 - P s + p 9 - P 10 

+Pn P 12 P 13 d~ P14 Pl5 T P 16 


r i0 = -Pi +P 2 +P 3 ~P4-Ps+Pe+P7 -Ps+P 9 - P 10 

Pll + P 12 P 13 P 14 P 15 P 16 

^ = -Pi +P2+P3-P4+P 5 -P 6 -P 7 +P 8 -P 9 + PlO 
TPn p 12 T P 13 P 14 P 15 T P 16 

r i2 = -Pi +p 2 -p 3 +p 4 +p 5 -p 6 +p 7 -p 8 -p 9 + p 10 
Pll P 12 P 13 P 14 d“ P 15 P 1G 

r i3 = -P 1 + P 2 + P 3 - P4 - Pg + Pg + P 7 - Pg - P 9 + PlO 

0-Pll P12 PlO + P14 + P 15 P 16 

r i4 = -Pi +P 2 -P 3 +P4-Po+P6-P 7 +P a +P 9 - PlO 
0-Pn Pi 2 T P 13 P 14 + P 15 Pig 


r =p —p —p +p — p +p +p —p —p +p 

15 *1 *2 l 3 *4 " 5 *6 l 7 l 8 *9 *10 

+p — p + p — p — p +p 

We interest to the positive partial transposition ( PPT ) 
region. If we consider the following notation 


Set ijP i ,P v P k iP l ) 


Pi +Pj+P k -Pi> 0 
Pi+P j ~P k +P l > 0 
Pi ~ Pi + P k + Pi > 0 
-Pi + Pi + P k + P, > 0 


then the positivity conditions for the eigenvalues of p TA 
are 

{ Set(p 1 ,p 2 ,p 15 ,p 16 ) 

Set (p 3 ,p 4 ,Pi 3 ,Pi 4 ) 

Set (p 3 ,p 6 , Pll ,p 12 ) 

Set (p 7 ,p 8 ,p 9 ,p w ) 

for the eigenvalues of p TB 

{ Set { Pl ,p 2 ,p 9 ,p 10 ) 

Set (p 3 ,p 4 ,p 11 ,p 12 ) 

Set (p 5 , Pel Pi 3 ,Pi 4 ) 

Set {p 7 ,p 8 ,p 15 ,Pi 6 ) 


( 10 ) 
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for the eigenvalues of p Tc 

{ Set (p 1 ,p 21 p 5 ,p B ) 
Set (p 3 ,p 4 ,p 7 ,p s ) 
Set (p 9 ,p 10 ,p 13 ,p li ) 
Set (jDjj , p 12 , p 15 , p lg ) 

for the eigenvalues of p T ° 

{ ^et(p 11 p 2 ,p 3 ,p 4 ) 
Set (p 3 ,p 6 ,p 7 ,p s ) 
Set^g,^,^,,^,) 
Set (p 13 , Pli ,p 15 ,p 16 ) 

for the eigenvalues of p TAB 

{ Set {p 1 ,p 2 ,p T ,p s ) 
Set (p 3 ,p 4 ,p 5 ,p B ) 
Set (p g , p 10 ,p 15 ,p 16 ) 

Set 


( 11 ) 


( 12 ) 



FIG. 1. The PPT feasible region for four qubit Bell diagonal 
states. For all points in the shadow region all eigenvalues of 
the density matrix ©> and all partial transposes are positive. 
The region is for inequalities m and urn 


for the eigenvalues of p TAC 

{ Set ( Pl ,p 2 , Pll ,p 12 ) 
Set (p 3 ,p i ,p g ,p w ) 
Set(p 5 ,p 6 ,p 15 ,p 16 ) 
Set(p 7 ,p 8 ,p 13 ,p 14 ) 

and finally for the eigenvalues of p TAD 


Adding both sides of above inequalities together and not¬ 
ing + ... + p 16 = 1, yields the following inequality 


(14) 

8 Pi - 6p 2 < 1 

(16) 

similarly, 

8 P 2 ~ bp, < 1 

(17) 


Set(p 1 ,p 2 ,p 13 ,p li ) 
Set(p 3 ,p 4 ,p 15 ,p 16 ) . . 

Set (p 5 ,p 6 ,p 9 ,p 10 ) 

Set(p 7 ,p 8 ,p 11 ,p 12 ) 

The total set of the above 4x4x7 = 112, inequalities 
define the PPT region of the four qubit Bell diagonal 
states ©• 


This region is illustrated in Fig[l] According to Fig{TJ if 
we consider vertex points (\,\) , (|,0) ,(0, |) and (0,0), 
because this points satisfy PPT conditions and feasible 
region is convex, then all of the points inside the shape 
fulfill the PPT conditions. 

We can find other feasible regions in the {p 1 ,p 3 ) plane, 
concerning the following inequalities 


A. Feasible regions 

Let us consider the feasible region which satisfy PPT 
conditions with details. To do this we can consider pj’s , 
that was ordered into eight pairs 


(Psi-PJ: (P 5 ,P 6 ), (Pr>Ps)> 


' Pi < P 2 + Ps+ Ps 
Pi < P 2 + P 7 + Ps 
< Pi < P 2 + P 9 + Pio 
P 3 < P 4 +Pn +P12 
P 3 < P4+P13+P14 
. P 3 < Pi +Pl 5 +PlB 

similarly we obtain 


=> 4( Pl +p 3 ) < 2 {p 2 +p 4 ) + 1 


(18) 


(PgjPio)’ (Pll’Pl 2 )l (Pl3’Pl 4 )> (P 15 ’ Pl6) 

Now if we choose one pair such as (p x ,p 2 ), then we can 
specify the feasible region in (p x , p 2 ) plane with the fol¬ 
lowing three inequalities ( PPT conditions) 

' Pi < P 2 + P 3 + P 4 
Pi < P 2 + P 5 + P 6 
Pi < P 2 + P 7 + P» 

< Pi < P 2 + P 9 + Pio 
Pi < P 2 + Pn + Pi 2 
Pi < P 2 + P 13 + Pi 4 
,Pi < P 2 + Pis + Pis 


4 (P 2 +P 4 ) < 2 (Pi +P 3 ) + 1 ( 19 ) 

from inequalities m and m we have 

4 (Pi +p 3 ) < 2{p 2 + Pi ) + 1 < Pl +p 3 + ^ (20) 

or 

Pi+P 3 <\ (21) 

so we have presented a new perspective from the spatial 
shape in two-dimension. 
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1. The first case 


2. The second case 


Let us consider the (yq , yq) and (yq,y> 4 ) planes. The 
PPT conditions for (yq,yq,yq,p 4 ) inequalities according 
to the following lines are satisfied 

/ 8p 1 - 6 yq = 1 f yq = 3yq/4 + 1/8 , , 

\ 8 p 3 - 6 yq = 1 \p 3 = 3yq/4 +1/8 

Similarly the PPT conditions for following cases 

(p 1 ,p 2 , P 5 , P e ) , (p l ,p 2 , P-n Ps) > (Pi,P 2 > P 9 ,Pia), 

(id 5 P 2 > Pll 5 Pl2 ) 1 (Pi 5 P‘2 1 PlS 5 Pl4 ) 5 (Pi 5 P2 > PlS > Pl6 )’ 

(jdsiPn P5 ’ Pe ) > (jdsiPn Pti Ps') i [PsiPai P 9 i Pio)> 

(.Psi P41P111P12) J (P3 i P 4 > F13 : P14 ) 1 (P3 : P 4 ; P15 i Pi6 ) 


are satisfied unless for the following cases 


’1 +P2 

+ Ps +Pe 

> 

0 

'1 +P2 

+ P 7 +Ps 

> 

0 

’1 +P2 

+ P 9 +P10 

> 

0 

’1 +P2 

+ P11 +P12 

> 

0 

’1 +P2 

+ Pis +Pia 

> 

0 

’1 +P2 

+ Pis +Pie 

> 

0 

’3 +Pa 

+ Ps+Pe 

> 

0 : 

’3 +Pa 

+ P 7 +Ps 

> 

0 

’3 +P4 

+ P 9 + P10 

> 

0 

’3 + P4 

+P11+P12 

> 

0 

’3 +Pa 

+ PlS + P 14 

> 

0 

’3 + P4 

+ PlS +Pl 6 

> 

0 

Ps +Pe +P2/ 4 

> 

1/8 

P 7 +P 8 + P 2 / 4 

> 

1/8 

P 9 +P10 + P 2 / 4 

> 

1/8 

P11 + 

Pl 2 +!+/ 4 

> 

1/8 

Pis + 

P 14 +P 2 / 4 

> 

1/8 

Pis + 

Pie + P 2 / 4 

> 

1/8 

Ps +Pe + P 4 / 4 

> 

1/8 

P 7 +Ps + P 4 / 4 

> 

1/8 

P 9 +P10 + P 4 / 4 

> 

1/8 

P11 + 

P12 + P 4 / 4 

> 

1/8 

Pis + 

P 14 + P 4 / 4 

> 

1/8 

Pis + 

Pie + P 4 / 4 

> 

1/8 


(23) 


(24) 


Adding the above inequalities yield to 

16 3 3 

2 Yl Vi + 2 ^ 2 +p 4 ) - 2 ( 25 ) 

i =5 

and using the normalization condition for yield 


2[ 1 -(Pi+P2+P3+P4)] + ^(P 2 +P 4 ) > ^ 


-^(P2+P 4 ) > 0 
(26) 


so we must have p 2 = p 4 = 0 =>• p x = p 3 = 4 . So the 
PPT conditions take the following simpler form 


' Ps +Pe > 1/8 
P 7 +Ps > 1/8 
P 9 +P 10 > 1/8 
Pn + P 12 > 1/8 
Pi 3 +P« >1/8 

. Pis + Pie > 1/8 


(27) 


In this case we consider if p 2 =0 and 8 p 3 — 6 p 4 = 1, 
then 

Pi ~Ps +Pa > 0 =>Pi + ^ > £ (28) 

so the PPT conditions for (yq, p 2 ; P 3 ■, Pa ) are satisfied. 
Adding the following PPT conditions 


-Ps 

+P4 

+ Ps 

+ Pe 

> 

0 

r Ps 

+ Pe 

+ P4/ 4 

> 

1/8 

-Ps 

+P4 

+p 7 

+ Ps 

> 

0 

p 7 

+ Ps 

+ F4/ 4 

> 

1/8 

-Ps 

+P4 

+p 9 

+ P10 

> 

0 

p 9 

+ P10 

+ P 4 / 4 

> 

1/8 

-Ps 

+ P4 

+ P11 

+ P12 

> 

0 ^ 

Pn 

+ P12 

+ F 4 / 4 

> 

1/8 

-Ps 

+ P4 

+ Pis 

+ P14 

> 

0 

Pis 

+ P14 

+ P 4 / 4 

> 

1/8 

-Ps 

+ P4 

+ Pis 

+ Pie 

> 

0 

- Pis 

+ Pie 

+ P 4 / 4 

> 

1/8 


(29) 


gives 

Ps +-+Pie + ^Pi > \ ( 80 ) 

Now if we set p x = | — ^ +e, where e > 0, the normaliza¬ 
tion condition, Yli=iPi = 1 ; and P 2 = 0 , 8 p 3 — 6 p 4 = 1 
yield to 

3 3 

^ + ^P 4 +Ps+Pe+ -+Pie = 4 ( 31 ) 

by applying (l30ll . we have —e>0=>e = 0, and p 4 = 
i-^-. Finally 


' Ps 

+ 

Pe 

+ P4/ 4 

= 1/8 

P 7 

+ 

Ps 

+ yq/ 4 

= 1/8 

P 9 

+ 

Pie 

+ P4/ 4 

= 1/8 

Pn 

+ 

P12 

+ F4/ 4 

= 1/8 ^ 

Pis 

+ 

P 14 

+ Jh/ 4 

= 1/8 

. Pis 

+ 

Pie 

+ F4/ 4 

= 1/8 

Pi 

— 

Ps 

+ Pe = 

= P 7 +Ps 


= 

P 9 

+ P10 = 

= Pll +P12 


= 

Pis 

+ P14 = 

= Pis +Pie 


(32) 


(33) 


B. MUB diagonal states which the PPT criterion 
is necessary and sufficient condition for separability 

In this section we investigate some MUB diagonal 
states which the PPT criterion is the necessary and suffi¬ 
cient condition for the separability of them. To this end, 
we write yq in the following pairs 


(PnP 2 ), (P 3 >P 4 )> (P 55 P 6 )> (P 7 >Ps)> ( 34 ) 

(PsiP io)> (Pn iPi 2 )i (^ 135 ^ 14)5 (Pib'iPie) 

Note that when any pair is zero then the two components 
of others are equal and the PPT criterion is necessary 
and sufficient for separability of the MUB diagonal den¬ 
sity matrix. For example, if we set yq = p 2 = 0 in the 
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first pair, then p 3 = p 41 p 5 = p 6 ,p 7 = p a ,p 9 = p 10 ,p n = 
P121P13 = Pu an d p 15 = p 16 • Now we can write the MUB 
diagonal density matrix in the following form 


1 

P= o 


2 iP3 (l^s) (V^3 I + I ^4 ) ("04 I ) + 2 ^ 5 (I "05 ) (^5 I + I ^6 ) (-06 I ) 

+ 2^7 (I ^7 ) (V’t I + I V’s } (V’s I ) + 2p g (I V’g } (V’a I + I V'lO ) (V’lO I ) 

+ 2 P 11 (| V’ll ) (Vhl I + 1 ^ 12 X ^12 I ) + 2 Pl 3 (I ^13 ) (^13 I 

(35) 


+ 1^14) (^14 I) + 2 P 1B (|^ 15 )(^15 I + IV’leXV’le I) 

which is a separable state. 


the Tr(Wp s ), takes the following simpler form 
Tr(Wp s ) = Aq ± Bo cos 9, cos 9 2 + 
sin 0, sin 9. sin 9, sin 9,. 


cos(^ 3 - p 4 ) 

( A\ cos p, cos p 2 + A 2 sin p, sin p 2 ) + 
sin(^ 3 + p 4 ) 

(A3 cos p, sin p 2 + A4 cos p 2 sin p, ) | (38) 
If we define new parameters 


III. WITNESSES DETECTING BOUND MUB 
DIAGONAL DENSITY MATRICES 


K 

h 3 


A + A 2 
2 

A + A 

2 


h 2 

K 


A ~ A 

2 

4 -A 

2 


(39) 


We introduce our linear four qubit entanglement wit¬ 
nesses that have the following generic form 


W = A 0 IIII ± B 0 a z a z II + 

A, (o’x&x&x&x “i” Gy ) -f- 

A 2 (<TyCTy<7 X <J X + OyOyOyOy) + 

A 3 (<J X (Jy (Ty U x + G x (JyG X G y) + 

A 4 ((Jy(J X Uy(J X - | Gy CTj G X Gy ) ( 3 G ) 

In order to investigate whether such operator really is an 
entanglement witness we must first prove its expectation 
value over separable states is nonnegative. To do so we 
evaluate the trace of witness over a pure product state 
[? ] which for four qubits state may be written as p s = 
|a)(a| ® |/3)( f3\ ® |y)(yI ® |A)( A|. The trace takes the 
following form 


Tr(Wp s ) = A 0 ± B 0 a 3 b 3 + 

A ( a 1 b 1 c 1 d 1 + a 1 b 1 c 2 d 2 ) + 

A 2 ( a 2 b 2 c,d 1 + a 2 b 2 c 2 d 2 ) + 
A 3 {a 1 b 2 c 2 d 1 + a 1 b 2 c 1 d 2 ) + 
A 4 (a 2 b 1 c 2 d 1 +a 2 b 1 c 1 d 2 ) (37) 


where 


Tr(\a)(a\(Ji) = a.i, Tr(|/3)(/?|cr, ; ) = 6 i; 
Tr(\i)('y\<Ji) = Ci, Tr(|A)(A|cr.j) = d. 


then 

Tr(Wp s 


= Aq± Bq cos 9, cos 9 2 + 
sin 9, sin 9 2 sin 9 3 sin 9 4 
cos(^ 3 - p 4 ) 

K cosip, -p 2 ) + h 2 cosip, + <P 3 ) 
sin (</> 3 + p 4 ) 

h 3 sin(y) 1 + p 2 ) + h 4 sin(< / o 2 - p,) 


By appropriate choice of the angles, one can minimize 
above expression, where its minimum value must be zero. 
For this purpose we set 9 3 = 9 4 = p 3 = p 4 = j , and 
define new parameters 


cos V’i 


sini/»i 


hi 


hi 


vw 

/14 

\/h\ + hi 


cos 


sin^ 


h 2 

h3 

\AiTh| 


then 

Tr(Wp s ) = = A 0 ± Bo cos 9 , cos 9 2 + sin^ sin 0 2 1 
\Jh\ + hi cos(V’i - pi + p 2 ) + 

\jh\ + h\ cos(V’2 - Pi - </>2)} 


Setting ipi = pi — p 2 and ip 2 = pi + P 2 , 


for * = 1,2,3 and cq’s are spin 1/2 Pauli matrices. With 
definitions 

a, = sin9 1 cosp 1 , a 2 = sin9 1 sinp 1 , a 3 = cos9, 

b, = sin9 2 cosp 2 , b 2 = sin9 2 sinp 2 , b 3 = cos9 2 

c, = sin9 3 cosp 3 , c 2 = sin9 3 sinp 3 , c 3 = cos9 3 

d, = sin9 4 cosp 4 , d 2 = sin9 4 sinp 4 , d 3 = cos9 4 


Tr{Wp s ) = A 0 ± B 0 cos 9, cos 9 2 + 

sin 9 1 sin 9 2 1 h\ + hi + \Jh\-\-h 3 1 
Using the identity 

— \Jr \ 1 + S 2 < ij cos 9 2 + S sin 9 2 < \/rj 1 + S 2 
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where p and 5 are coefficients of cos O 2 and sin 62 respec¬ 
tively, we have 

Tr{W Ps ) > A o = F {b o 2 cos 2 0 1 + 



choosing 

Ao = B 0 = jh* + hi + yjhl + hi (40) 

yields to 


Tr{Wp s ) > 0 

and the entanglement witness becomes 
r a 

W = A 0 IIII ±a z a z II+— L (a x a x a x a x +a x a x a y a y )+ 

L A 0 

A A 

— ((JyCTyCT X (T X A Cf y(J y(7y<7y) -\ “ [(7 X (7y(T y(7 X U O X CI y<J X (Ty)-\- 

Aq A 0 

A 4 , \ 

\(JyO'x&yGx\&y&x&x&y) 

^0 


we note that 


A 1 _ h x + h 2 _ cos i/j 1 y/h% + h'l + cos ip 2 y/h\ + 

Aq Aq Aq 

and if 

p =i\ / ^=2i\/ 2< ^ + - 4 ‘ ) (41) 

then we have 



FIG. 2. Non linear EW is the envelope of a family of lin¬ 
ear EW s. The gray region is the PPT entangled or feasible 
region. Each linear EW separates the separable region from 
the entangled one. The set of such linear EW s is a family of 
curves which their envelope can be considered as a nonlinear 
EW. 


construct nonlinear entanglement witnesses for the four 
qubit MUB diagonal states using the envelope of a family 
of curves definition. 

Let F : 5ft x Sft r —► 5ft be a smooth map and t. cci, x r 
coordinates on the left. Consider F as a family of func¬ 
tions of x, parameterized by t. The envelope, of the fam¬ 
ily F is the set |38| 

PLf = { 2 ; £ 5ft’ : there exixts t £ 5ft with (43) 
F(t, x ) = dF{t, x)/dt = 0} 

Using this definition we can find the envelope of our linear 
EWs. This envelope corresponds to a nonlinear EW. To 
do so, we consider the trace of EW over four qubit MUB 
diagonal state, ©, as a family of functions (linear EWs) 


—A = p COS7/q + (l— p) COS 1p 2 , -A- = pCOS^q —(1—p) COS^q 

Aq a 0 


~T =P sinV> 1 +(l-p)sinV> 2 , -p = -psin^ 1 +(l-p)sinV > 2 

Aq A 0 

and the entanglement witness, W, can be written as the 
following form ( without loss of generality we divide the 
W by A 0 ) 

W = IIII ±<r z a z II + 

[pcos Vh + (1 - p) cos tp 2 j (a x <j x a x cr x + a x a x a y a v ) + 
[pcosi/q — (1 — p) cost/> 2 ] ( a y a y a x a x + a y a y a y a y ) + 
[p sin ip 1 + (1 - p) sin ijj 2 ] [a x a y a y a x + a x a y a x a v ) + 
[- pshn/q + (1 -p) sin i [) 2 ] (a y a x a y a x + a y a x a x a y ) 

(42) 


Tr(Wp) = 1 ± r g + (jp cos V*! + (1 — p) cos?/q,)(r 8 + rqj 
T(p cos ip 1 - (1 -p) cos ij) 2 ) (r 14 +r 1B ) 

+ (psinV> 1 + (1 ~p)smi/j 2 )(r 9 +r 10 ) 

+ (-p sin ip 1 + (1 - p) sin tp 2 ) (r 12 + r 13 ) (44) 


or in terms of pi 


16 


Tr(Wp) = 1 ± [ l-2^p, 
3 =9 


4p 


(p — p + p — p )cost/>i 

V/ ll ^ 12 “ 13 *14 ' ^ 


+(P 9 -Pio~P 15 +Pj sin ^ 


+ 


4(1-p) (p -p +p -pJcosV >2 


- (p i - P 2 - P 7 + P s ) sin i >2 


(45) 


IV. NON-LINEAR ENTANGLEMENT 
WITNESSES 


This family of functions has two parameters, ?/q and ip 2 
and the condition dTr(Wp)/dip 1 = 0 yields 


We showed that how one can find the PPT feasible 
region and introduced the linear EWs. Now we can 


ip 1 = arctan 


( P 9 -P 10 ~P 15 +P, 6 \ 
Vn-Pu + Pis-PiJ 


(46) 




















similarly dTr(W p)/dip 2 = 0 leads to 

^~P 1 + P 2 + P 7 ~P 8 


0 2 = arctan 


p-p+p-p a 


(47) 


now if we insert equations (l46l) and (l47l) in (l44l) and sim¬ 
plify the result then 


16 


Tr(Wp) = l + a 0 (l-2j> 

i=9 


+ 4 P (p 41 - p 12 + p 13 - p 14 ) 2 

1 / 2 

4”(Pg Pic P 15 4” Pie) 
x sgn (p 14 -p 12 + p 13 -p 14 ) 

+ 4(l-p) a 2 (p 3 -P 4 +P 5 -P 6 ) 2 

-I 1/2 

+ (Pi -P 2 -P 7 +Ps) 2 
x sgn (p 3 -p 4 +p 5 -p 6 ) 


eg 

a 


(48) 



where a 0 = ±1, sgn is the sign function, p = 

—!j- y/ 2 (A^ + A^), and if 2 nk — f < 0 4 < 27rfc + ^ then 
ai is +1, otherwise aj, is —1 for i = 1,2. Now the enve¬ 
lope equation is Tr(Wp) = 0. The following example is 
given to indicate the full idea and details of the envelope 
approach for nonlinear witnesses. 


FIG. 3. Linear Entanglement Witnesses. Any point in the 
shaded region represents a possible four qubits density ma¬ 
trix. Each line shows a linear witness separating the entan¬ 
gled density matrices from separable. Here we plotted 50 
linear witnesses for various values of § < 02 < yy. 


Example 


Let us consider a simple example of envelope algorithm 
in operation for a special case where 


p = 1/10, p i = 0 for i > 4, p 3 = 1 — p 3 — p 2 

From equation (ll5l) and taking the plus sign in the second 
term, we have 


18 18 

Tr(Wp) = 2+y(p 2 -pi)sinV>2 + y(l-pi—P 2 )cosi /2 


(49) 

where P 2 , P 2 are variables of the density matrix with 
constraints 0 < p\ < 1, 0 < P 2 < 1, 0 < p\ +P 2 < 1, and 
0 < 02 < 27 t, is a witness parameter. Solving Tr(Wp) = 
0 in terms of P 2 yields to 


P2 


cos 0 2 - Pi (cos ip 2 + sin 02 ) + 5/9 
cos 02 — sin 02 


(50) 


Now any value of 02 corresponds to a linear witness. 
Fig. [3] shows 50 of these linear witnesses for 02 = 
{f’ f + I)’ f + !§’■• • i if }’ in the region restricted to 
constraints 0 < pi < 1, 0 < p 2 < 1, and 0 < p\ +p 2 < 1. 
We can find the envelope equation of these linear wit¬ 
nesses using definition, (l43l) . If we obtain 0 2 from 
dTr(Wp)/di/j 2 = 0 and insert it in Eq. (IfiUl) then 


P2 = ^ (9 - ^324^(1-pi)-3l) 


for the interval, ^ (9 — 5\/2) < pi < | , and 


P2 = ^ (9 + y/32A Pl {l - pi) - 3l) 


for the interval, yg (9 — 5-02) < pi < |. 

This is the envelope equation with constraints 0 < 
Pi + P 2 < 1. As this is a nonlinear function we call 
such envelope as nonlinear entanglement witness. The 
PPT region, the detected entangled region, and the non¬ 
linear EW as the envelope are plotted in Fig. [4] The 
PPT region has obtained from (flTH) and 03 which is 
8 p 1 — 6 p 2 < 1 and 8 p 2 — 6 p 1 < 1 ( See Fig. |T] ). As can 
be seen the nonlinear witness has wider detection range 
of entanglement. The region enclosed with the bold lines 
corresponds to the positive partial transpose entangled 
states (density matrices) which can be detected by non¬ 
linear witnesses, that is to say, Tr(Wp) < 0 . 

Consequently, we achieved the nonlinear EWs for some 
four qubit MUB diagonal density matrices and the neg¬ 
ativity of (H5|) is the evidence of four qubit entanglement 
in the system. At the end we present other nonlinear 
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Pi 


V. THERMAL ENTANGLEMENT AND STATE 
PREPARATION FIDELITY 


Here we consider a canonical ensemble of four qubits 
identical systems in thermal equilibrium and we find the 
entanglement detection condition at temperature, T. We 
will show while the fidelity criterion fails to detect the 
entanglement, our nonlinear witnesses can detect the en¬ 
tanglement. 

From (l48l) one can find the entanglement dependence 
on the temperature. Rewriting this condition in terms of 
the temperature (here we set ks = 1), 


Tr(Wp ) = 1 + a 0 (l - | e ~ Ei/T ) + 
\ z j=g J 


Apa 1 


l ~ E n 
[e t 


T — p T 


-Eq -E I 


e T - e T - e T -f e T 


1/2 


FIG. 4. The PPT region (the dotted light gray region, see also 
Fig. m. detected entanglement ( the dark gray region), and 
nonlinear EW as the envelope ( the dashed curve ) is plotted 
for special case which a 0 = +1, a 1 = +1, a 2 = —1, p = 
Tj , p. = 0 for i > 4. The nonlinearity of EW is obvious and 
in this case is a semicircle. We see that the nonlinear EW 
detects more entangled density matrices. 


EW s families. These have the following form 
W {il , i 2 ,i 3M} = IIII±O j 


+ 

Ox&x@x&x ( 1) 1 ®x 

Ai/Aq 

+ 

®x®x ( 1) 

A 2 /A 0 

+ 

&x&y&y@x “1“ ( 1) 3 &x&y&x&y 

Ag/A 0 

+ 

@y®x®y®x H - ( 1) 4 GyGx®x®y 

Aa/Ad 


(51) 


here 


{*1, *2j *3, h} G 

{{ 0 , 0 , 0 , 0 }, { 0 , 0 , 1 ,!}, { 1 , 1 , 0 , 0 }, { 1 , 1 , 1 ,!}}, 


and 


Oj G a z a z II,IIa z a z ,Ia z Ia z ,Ia z a z I , 

(7 Z (J Z II, CT Z I<T Z I, (T Z II(J Z , (T Z (T Z <J Z (7 Z 


Similarly one can show that all 64 relations of m are 
nonlinear EW s. Furthermore if we consider the notation 
P(m, n) for permutation of mth and nth Pauli matri¬ 
ces of the eight terms in the brackets of then the 
permutation P(l, 2), P(l, 3), P(l, 4), P(2, 3), and 
P(3, 4) gives new nonlinear EW s. Therefore we have in 
total 64 x 5 = 320 nonlinear EWs. 


x sgn 


~ E 11 ~ E 12 — e 13 ~ e 14 

' e t — e T + e t — e t 


4(1 — p)a 2 


^3 — E 4 —-^5 _ 

e t _ e T + e T — e T 


( zE 1 

(e * 


— e 2 —E t — \ 2 

— e t — e T + e T 


x sgn( 

in this expression 


— e 3 —E 4 ~ e 5 ~ e 6 

' e t — e t + e t — e t 


Z 


16 

z = J2 e ~ ¥ 

3 = 1 


+ 


(52) 


(53) 


is the partition function of the system. For a given Ej s, 
the negativity of (1521) for a temperature interval, is the 
sufficient entanglement condition. 

For example we consider the following Hamiltonian 

H = — J (a z a z II + Ia z a z I + IIcr z a z 

+a z IIa z + cr z Ia z I + Ia z Ia z ^j 

~\~h^(T x (J x (J x (Jx T T &z&z&z&z'j ( 54 ) 

the hrst part shows the pair coupling between four 
qubits with coupling constant J, the second part rep¬ 
resents some spin interactions among all four qubits with 
strength h. This Hamiltonian in the Bell-diagonal bases 
is diagonal with energy eigenvalues 

Ei = 3h — 6J, E 2 = —h — 6J, 

Eg = E4 = E 5 = Eq = Eg = E 10 = P15 = Ei e = —h, 

E’j = En = Pi 3 = 3/i + 2 J, 

Eg = E 12 = E 14 = —h + 2 J 


( 55 ) 
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FIG. 5. The entanglement detection in terms of temperature, 
T for three values of J = — 1, —2, —3, —oo. Entanglement is 
detected for Tr(Wp) < 0. Here we choose ao = a\ = p = h = 

1. 



FIG. 6. The fidelity as a function of temperature for Hamil¬ 
tonian (15411 . with h = 1. Please note F(T) —> 0 as J —> oo. 
There is no entanglement detection at any temperature. 


TABLE I. Threshold temperature for Hamiltonian (1541) . 
J -1 -2 -3 -oo 

T th 2.6135 3.6232 4.3418 7.8305 


and the partition function of the canonical ensemble of 
four qubits systems at equilibrium temperature, T, is 


h — 3h — 2 J 

Z = 8c T + 3c 


fa+6J . . 

e~ (56) 


Now one can use (l52l) for calculating the entanglement 
detection condition at temperature T. The Tr(Wp) is 
equal with 


1 + 


a 0 


(l +eT^j ^ 


-1 


8 J 

■ e T 


-4 h _ 2(2 h+J) 

3 + 3eT" + 8e— t — 


+ 8aip ^1 — e 4 r ^ - 

4h \ 8J 

[1 + eT jer 


(57) 


Using the negativity of this, we can estimate the thresh¬ 
old temperature, T t h, for entanglement in the thermody¬ 
namic limit. Fig. [5] shows the plot of this condition in 
terms of temperature for ao = a\ = p = h = 1, and three 
values of J = —1, —2, —3, —oo. For T < T t h , the entan¬ 
glement has been detected by our nonlinear witness (1481) . 
These threshold temperature is represented in Table. nit 
is seen that threshold temperature increases and reaches 
to its maximum with the decrease of J. 

It is interesting to consider the detection condition for 
extreme limits of J and h, 


lim Tr(Wp ) = 

h—>oo,J—} — oo 


1 - -(ao + 8aip) 


< 0 


(58) 


or (a o + 8 aip) > 3. This result is independent of temper¬ 
ature and for values such as ao = a\ = p = 1, not only 
the system is entangled in any temperature but also we 
have a witness detecting this entanglement. 


Now we compare the state preparation fidelity F of a 
iV-qubit state to investigate the four qubits Bell diagonal 
state entanglement. The fidelity is defined as [39] 

F(p) = (i/; 

GHZ I/# GHZ ) (59) 

where for four qubits, 1 4> GHZ ) = (|0000) + 11111 ))/a/2 , 
and p is our four qubit Bell diagonal density matrix. The 
sufficient condition for four qubits entanglement is phi . 

S3 


F(p) > i (60) 

Using definition (l59ll we calculate the fidelity for the pre¬ 
vious example with h = 1. The result is 


F(T) = 


e t 


o . O 4 2(2 +J) 8J 

3 + 3er+8e t +er 


^1 + er j 


(61) 


Fig. \E\ shows the fidelity for some values of J. Clearly 
fidelity is zero as J —> — oo. We see that for negative val¬ 
ues of J , limU(T) —>1/16 when T —> oo. In other words, 
the fidelity criterion fails to detect the entanglement for 
this example. 


VI. CONCLUSION 

We have constructed linear and nonlinear entangle¬ 
ment witnesses with wider detection region for four 
qubits systems in mutually unbiased bases. These wit¬ 
nesses may detect the entangled positive partial trans¬ 
pose density matrices. We have extended them to a 
canonical ensemble consisting of infinite number of four 
qubits in thermal equilibrium and have detected the en¬ 
tanglement at any temperature for a special case. Results 
shows that these witnesses also can detect the entangle¬ 
ment in an ensemble consisting four qubits systems for 
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some Hamiltonians while another method, state prepa¬ 
ration fidelity, fails to detect. 

Our results highlight the potential of this method and 
may be of considerable practical value. Building on this 
success one can begin other trial systems for entangle¬ 
ment investigation. Accurately observing entanglement 
is a demanding task from the experimental point of view 
and further studies will probe fundamental aspects of 
usefulness of nonlinear EWs. In this area of study much 


practical research such as detecting entanglement in tele¬ 
portation, quantum cryptography, and quantum algo¬ 
rithms remains to be done. In conclusion, this paper 
has explored the nonlinear EWs family for detecting the 
PPT states of a given four qubit density matrix. Our 
study focus on the aspects of wider range detection of en¬ 
tanglement with nonlinear EWs and building them with 
linear EWs. 
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